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During the last few years there has been a renaissance in charm and charmonium spectroscopy 
with higher precision measurements at the ip and '0(3770) coming from BESH and CLEOc and many 
new discoveries coming from B-factories. In this paper, I review some new results on "classical" 
charmonium and e^e~ hadrons using B-factory Initial State Radiation and two photon events. 



I. INTRODUCTION 

\ In recent years, there has been tremendous progress in charm and charmonium spectroscopy with many new 
results from Belle, BaBar, CDF, DO, BES and CLEOc. I will review recent results on "classical" charmonium; 
O ■ Galina Pakhlova (see her paper in these proceedings) covers exotic charmonium states: X, Y, Z, etc. 

, Charmonium provides detailed information on QCD in the perturbative and non-perturbative regimes, as well as 
f — I providing a laboratory for precision tests of lattice QCD and effective field theory [1]. 

Charmonium may be produced by e+e~ annihilation, B decays, e+e~ annihilation where one or both electrons 
loses energy by Initial State Radiation (ISR), and two photon processes. The very high luminosity of B-factories 
^ \ allow the use of the latter two even though their cross sections are suppressed. Only J^^ = 1~~ charmonium states 
can be produced directly in e+e~ annihilation. However states below the ^1^(28) may be produced by ^p{2S) radiative 
and hadronic transitions, a technique which has been used extensively by BES and recently by CLEOc. 
. I will also cover some light hadron spectroscopy results, including the y(2175) and (/)(1680), two photon results at 
Belle, and ISR results at BaBar, and I will report the status of BEPCII and BESIII. I apologize to everyone whose 
results I do not cover because of lack of time and space. Please see the references for details. 
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II. CHARMONIUM RESULTS 

Recently CLEOc obtained 27 million '0(25') events and has reported many new high precision results using this 
sample. This is the world's largest '0(25') sample produced in e+e~ collisions. 

A. i/?(2aS') — > J/'^ transitions 

CLEOc has studied il^{2S) J/ijj transitions using J jil) ^ e+e~ and \i~ decays to identify the J/'0 [2]. A 
summary of their branching fraction results is shown in Table I, along with those of the Particle Data Group (PDG). 
Using the product Xc branching fractions and B{ip{2S) jXcj) = (9-4 ± 0.4)%, (8.8 ± 0.04)%, and (8.3 =b 0.4)% [3] 
for J = 0, 1, 2, respectively, the B{xcJ 7^/0) radiative branching fractions are obtained and also listed in Table I. 
Precision branching fractions are important since 0(25) production and decay is a primary means of producing Xc 
states, and 0(25) production in many experiments is measured by the narrow mass peak recoiling against the tt+tt" 
in 0(25) TT+TT" J/0 decay. 



B. J/i/j and V^(25) -/rjc 

The precise determination of the rjc mass provides information on the hyperfine splitting of the rjc and J/0. 
However, although there have been many measurements of the r]c mass and width, the measurements do not agree 
very well, and the fitted masses and widths in the PDG [3] have very low confidence levels: 0.002 for the mass and 
< 0.0001 for the width. Further, the masses obtained from J/0 and 0(25) decays are about 5.3 MeV/c^ or 3cr 
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TABLE I: Branching fractions for '0(25') J/ip transitions from CLEOc [2]. B{iIj{2S) IXcj) values from Ref. [3] are used 
to obtain the CLEOc Xc lJ 1'^ radiative branching fractions, and the last error for them is from this branching fraction. 



Channel 


B(%) CLEOc 


B(%) PDG08 [3] 


Tl^'K~ J /l^ 

7Xco(XcO Jl"^) 
7Xci(Xci J/i^) 

7Xc2(Xc2 Jj"^) 

anything J/'0 


35.04 ±0.07 ±0.77 
17.69 ±0.08 ±0.53 
3.43 ± 0.04 ± 0.09 
0.133 ±0.008 ±0.003 
0.125 ± 0.007 ± 0.013 
3.56 ±0.03 ±0.12 
1.95 ±0.02 ±0.07 
62.54 ± 0.16 ± 1.55 


32.6 ±0.5 
16.84 ±0.33 
3.16 ±0.07 
0.126 ±0.013 
0.120 ± 0.010 
3.15 ±0.08 
1.66 ±0.04 
57.4 ± 0.9 


Xco iJ/i^ 
Xci iJ/i^ 
Xc2 iJ/i^ 


1.32 ±0.07 ±0.14 ±0.06 
40.5 ±0.3 ±1.4 ±1.8 
23.5 ±0.2 ±0.8 ±1.1 


1.28 ±0.11 
36.0 ± 1.9 
20.0 ± 1.0 



lower than those obtained from 77 fusion and pp annihilation. Another problem is that the branching fraction for 
J jil) ^ jr]c^ B{J /il) jrjc) = 1.3 ± 0.4% [3], is very low compared to recent Lattice QCD results, B{J/tlj jrjc) = 
2.1 ±0.1 ±0.4% [4]. 

Studying J/i/j ^ jrjc and i^{2S) jrjc is very important since these are magnetic dipole (Ml) transitions, and 
their branching fractions are necessary for normalizing rjc branching fractions. 

CLEOc studied these decays using 24.5 million '0(26') events [5]. They obtained the branching fractions by fitting 
the 7 energy spectrum using three samples, (1) '0(26) jrjc inclusive, (2) '0(26) tt+tt" J/'0, J/ip ^ ji^c^ Vc 
and (3) ^^(26) jrjc^ rjc X^, where Xi denotes 12 exclusive rjc hadronic decays. Sample (2) was used to determine 
the 7 line shape, shown in Fig. 1. Interestingly, they find that the line shape can not be fitted by a simple Breit- 
Wigner plus a resolution function, and an empirical shape was used. Their branching fraction results compared with 
the PDG are shown in Table IL The values are very different than the PDG and will affect all rfc branching fractions. 
Their B{J/ip jrjc) agrees much better with the Lattice QCD prediction above. 
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FIG. 1: Gamma energy spectrum from CLEOc [5] for J/ip 'jric using '0(26) tt^tt" J/'0, J/ip jric, rjc Xi events, 
where Xi denotes 12 exclusive rfc hadronic decays. A fit using a Breit-Wigner plus resolution function (dotted curve) is unable 
to fit the data (dots with error bars). 

If they fit the three samples with an unmodified Breit-Wigner, they obtain a mass m{r]c) = 2976.7 ± 0.6 MeV/c^ 
(stat. error), compared to m{r]c) = 2982.2 ±0.6 MeV/c^ with their empirical line shape. The first mass is consistent 
with other determinations from J/ip and il){2S) decay, while the second is consistent with those coming from 77 
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TABLE II: Branching fraction results for ip{2S) and J/ip ^ jrjc transitions from CLEOc [5]. 



Channel 


CLEOc 


PDG08 [3] 




(4.32 ±0.16 ±0.60) X 10"^ 
(1.98 ± 0.009 ± 0.30)% 


(3.0 ±0.5) X 10-^ 
(1.3 ±0.4)% 



fusion and pp annihilation. As examples of recent 27 production m{rjc) measurements. Belle determined m{r]c) = 
2986.1±1.0±2.5 MeV/c^ in r/c ^ four body decays [6] and m{r]c) = 2981.4±0.5±0.4 MeV/c^ using r]c KsKn [7]. 
The line shape problem may explain the 3a mass difference with 77 fusion and pp annihilation, but the uncertainty 
in how to deal with the line shape, according to the authors, prohibits the precise determinations of the rjc mass and 
width. The authors point out that understanding the energy dependence of the ip{lS, 2S) — jr]c matrix element is 
crucial for understanding radiative decays. 

Christine Davies showed a comparison of Lattice QCD (2007 HPQCD/MILC/FNAL) results with experiment [8]. 
One of the biggest discrepancies was in the m{J/ip) — m{r]c) comparison. Assuming from the CLEOc mass mea- 
surements with and without using the empirical line shape that the average r]c mass might shift upwards by roughly 
3 MeV/c^, the agreement between LQCD prediction and experiment for m(J/'0) — m{r]c) would be considerably 
improved. 

C. hcCPi) 

In 2005, E835 [9] and CLEO [10] reported measurements of the mass of the hd^Pi). CLEO used e+e" ij{2S) 
TT^hc and determined the he mass by measuring the mass recoiling against the tt^ using both he 'jrjc inclusive events 
and exclusive rjc decays. They have repeated their analysis using the 25 million il){2S) sample [11]. They find excellent 
agreement between the inclusive and exclusive results and obtain a combined result of m{hc) = 3525.28 ± 0.19 ± 0.12 
MeV/c^. Combining with their 2005 result, they obtain m{hc)AVG = 3525.2 ± 0.18 ± 0.12 MeV/c^ and a product 
branching fraction, {Bi x B2)avg = (4.16 ±0.30 ±0.37) x 10~^. A precise determination of the mass is important to 
learn about the hyperfine (spin-spin) interaction of the P wave states. Using the spin weighted centroid of the ^Pj 
states, < m{^Pj) >, to represent m{^Pj), they obtain Am/,/(lP) =< m{^Pj) > -m(^Pi) = ±0.08 ± 0.18 ± 0.12 
MeV/c^. This is consistent with the lowest order expectation of zero. 

D. XcJ 77 

Xcji^Pj) 77 decays (QCD) are analogous to triplet decays of positronium (QED). For R = T{^P2 
77)/r(^Po 11)^ even the differences due to different masses and wave functions cancel, and for both R = 4/15 ~ 
0.27 [12]. Departures from this are from strong radiative corrections and relativistic effects. 

CLEOc has studied this process using tp{2S) JiXcJiXcJ 7273 with their 25 million event sample [13]. The 
71 energy spectrum is shown in Fig. 2, where clear peaks corresponding to XcO and Xc2 are seen. Xci is forbidden 
by the Landau- Yang Theorem [14]. Results are listed in Table IIL Averaging the value of R with those from 
previous experiments, < R >= 0.20 ± 0.02 is obtained. The theoretical first order pQCD prediction for R is 
Rj.^ = (4/15) [1 - 1.76as] [15], and for as = 0.32, 

Rrh = 0.12. 

The disagreement with the experimental result confirms the inadequacy of the first order radiative corrections. 

Belle has also determined R using a number of exclusive decays of the Xc's produced in 77 Xc0,2 [6]. They 
measure F^^ x B{xc0,2 -^), where X is an exclusive decay, and divide by the known branching fraction to obtain 
F^^. The results are also listed in Table III, and their R values agree very well with < R> determined by CLEOc. 
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FIG. 2: Energy spectrum of 71 in ip{2S) 7iXcJ,XcJ 7273 from CLEOc [13]. 



TABLE III: XcJ ^ 77 widths and R = r(^P2 



77)/r(^Po ^ 77) from CLEOc [13] for XcJ ^ 77 and Belle [6] for 
77 ~^ XcJ, with XcJ ^ exclusive hadronic decays. The CLEOc result is the first result; their last error is the contribution 
from S(^(25)^7Xc) [3]. 



Channel 




r77(Xc2) 


R 


XcJ 77 


2.53 zb 0.37 zb 0.11 zb 0.24 keV 


0.60 zb 0.06 zb 0.03 zb 0.05 keV 


0.237 zb 0.043 zb 0.015 zb 0.031 


^KsKs 

2K2ti 
^4K 


2.53 zb 0.23 zb 0.40 keV 
1.84 zb 0.15 zb 0.27 keV 
2.07 zb 0.20 zb 0.40 keV 
2.88 zb 0.47 zb 0.53 keV 


0.46 zb 0.08 zb 0.09 keV 
0.40 zb 0.04 zb 0.07 keV 
0.44 zb 0.04 zb 0.16 keV 
0.62 zb 0.12 zb 0.12 keV 


0.18 zb 0.03 zb 0.04 
0.22 zb 0.03 zb 0.05 
0.21 zb 0.03 zb 0.09 
0.21 zb 0.06 zb 0.06 



E. Anomalous line shape of a{e^e hadrons) in the i/?(3770) energy region 

BESII has accumulated more than 30 pb~^ of data in the region of the -0(3770) from 3.650 to 3.872 GeV, while 
CLEOc has 818 pb~^ at the ?/^(3770). These samples have provided precision measurements of D meson decays using 
the very clean '0(3770) DD events, as well as vastly improved knowledge about the '0(3770). However, there has 
been a puzzle concerning the amount of non- DD decay of the '0(3770). The 0(3770 is just above threshold for DD 
production and is expected to decay into DD pairs with a branching fraction greater than 98%. Surprisingly, BES 
measured the branching fraction of '0(3770 decays to DD to be B{ij{3770) ^ DD = (85 ±3)% [3, 16, 17] and directly 
measured 5('0(377O) ^ non-DD = (13.4 zb 5.0 =b 3.6)% [18] and 5('0(377O) ^ non-DD = (15.1 =b 5.6 =b 1.8)% [19]. 
However, BES and CLEOc have searched for exclusive non-DD decays of the 0(3770), and the summed non-DD 
branching fractions measured by each of the collaborations are less than 2% [20, 21]. 

In the inclusive measurements, BES assumed a single resonance in the energy region between 3.7 and 3.872 GeV. 
To understand the discrepancy between the inclusive and exclusive measurements, BES has reanalyzed the fine R- 
scan {R = ^'^^^^J^f^drons ^ ^Yii^ region and finds that the fit to a single resonance is very poor and that allowing 
two non-interfering or two interfering resonances gives a much better fit, as seen in Fig. 3 [22]. The large non-DD 
inclusive branching fractions measured by BES may be due partially to the assumption of only one simple resonance 
in this region. This cross section anomaly must be confirmed, and this will be a high priority for BEPCH and BESHI 
(see below). 
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FIG. 3: Inclusive hadronic cross section in the region of the '0(3770) by BES versus center of mass (CM) energy [22]. (a) Fit 
with one resonance and (b) fit with two interfering resonances. 



III. r(2175) AND 0(1680) 



A structure at 2175 MeV in the 0/o(98O) mass was observed by BaBar in the ISR process e+e — > 7^^^0/o(98O); 
the mass and width are m(r(2175)) = 2175 ± 10 ± 15 MeV/c^ and r(r(2175)) = 58 ± 16 ± 20 MeV/c^ [23]. BaBar 
speculated that the y(2175) is the ss version of the y(4260) [24] since it is also a 1~~ and has somewhat similar 
decay properties [23]. 

BES searched for the r(2175) in J/V^ ^ 7^0/o(98O), 7? ^ 77, ^ i^+K", /o(980) ^ tt+tt- using 58 million J/i/j 
events and found a peak in the (/)/o(980) mass around 2175 MeV/c^ [25]. Fig. 4 shows the simultaneous fit to signal 
and sideband events with a Breit-Wigner to represent the signal and a third order polynomial for the background. 
The peak has a significance of about 5 a, and the mass and width obtained are m(F(2175)) = 2186 ± 10 ± 6 MeV/c^ 
and r(y(2175)) = 65 ±23 ±17 MeV/c^, in good agreement with BaBar. Fitting also the smaller ~ 2a peak at around 
2460 MeV/c^, also seen by BaBar, does not change the mass and width of the first peak. The product branching 



fraction is B{J/^ r]Y{2175)) • 5(^(2175) ^ 0/o(98O)) • 5(/o(980) 



-) = (3.23 ±0.75 ±0.73) x 10" 




(b) 



J3.0tj 



^■I.USJ iJ.lt! :i,^<ii ^.3Ei i?^l5 



FIG. 4: Fit of y(2175) by BES [25]. (a) Fit (solid curve) to the data (points with error bars) versus m(0/o(98O); the dashed 
curve indicates the background function, (b) Simultaneous fit to the sideband background events (points with error bars) 
with the same background function. The background normalizations for the two plots are constrained to be equal. 



Belle has also searched for the ^(2175) in e+e" ^ 7,^^0/o(98O) using 673 fb-^ of data at the T(4S') [26]. They 
find a peak in the e+e~ (/)/o(980) cross section and fit it with one Breit-Wigner interfering with a non-resident 
background (see Fig. 5 b), as in the BaBar analysis, and also with two Breit-Wigner functions. They also see a peak 
in e+e~ 1isr^^^^~ (^^^ ^^^S- 5 a) and fit it. For their final results, they take a simple average of their fits and 
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enlarge errors to cover the spread. They find m(Y'(2175)) = 2133t??5 MeV/c^ and r(y2175)) = I69tlf MeV/c^, 
where the Belle errors include both statistical and systematic errors. 

Belle also finds for the first time in e+e~ 1isr^^~^'^~ ^ clear ^(1680) (see Fig. 5 a), which was first seen by DM1 
25 years ago [27]. Belle determines m((/)(1680)) = 1687 ± 21 MeV/c^ and r((/)(1680)) = 212 ± 29 MeV/c^ BaBar 
also recently reported the (/)(1680) in e+e~ 7isr^^ 7isr^^^ [2^]- "^^^ masses and widths are summarized in 
Table IV along with the PDG [3] 0(1680) values. The Belle results are preliminary. 



TABLE IV: Summary of masses and widths of y(2175) and 0(1680). 



Experiment 


Channel 


Mass (MeV/c^) 


Width (MeV/c^) 


BaBar [23] 
BES [25] 
Belle [26] 


y(2175) ^ (/)/o(980) 
y(2175) ^ (/)/o(980) 
y(2175) ^ 07r+7r-,(/)/o(98O) 


2175 ±10 ±15 
2186 ± 10 ±6 
2133t??5 


58 ± 16 ± 20 
65 ±23 ±17 
1691^°/ 


Belle [26] 
BaBar [28] 
PDG [3] 


0(1680) 07r+7r- 
0(1680) ^ K*K and 0ry 
0(1680) 


1687 ±21 
1709 ± 20 ± 43 
1680 ± 20 


212 ±29 
322 ± 77 ± 160 
150 ± 50 



Belle finds a wider F (2175) and notes that the widths of the F(2175) and 0(1680) are rather similar, which suggests 
the possibility that the y(2175) may be an excited state. So what is the y(2175)? It could be a ss analogue of the 
y(4160), as suggested by BaBar; a ssg hybrid [29]; a 2^Di ss state [30]; a ssss tetraquark state [31]; a AA state [32]; 
or, as suggested, by Belle a conventional ss state. More data are needed to understand the 1^(2175). 




FIG. 5: Belle y(2175) and 0(1680) fits, (a) Fit to 07r+7r" cross section in e+e~ 7^5^07r+7r" and (b) fit to 0/o(98O) cross 
section in e^e~ 7^^^0/o(98O) with one Breit-Wigner (dashed curve) interfering with a non-resonant contribution (lower 
smooth curve) [26]. 



IV. BELLE 27 PHYSICS 

Two 7 collisions provide valuable information on both light and heavy quark resonances, perturbative and non- 
perturbative Quantum Chromodynamics, and hadron production mechanisms. Some 77 charmonium results have 
already been reported in Sections II B and II D. Recently Belle studied 27 production of tt+tt" [33], and now has 
new high statistics results on 27 production of tt^tt^ using 95 fb~^ of data [34]. 

Shown in Fig. 6 is a partial wave analysis fit of the differential cross section in the low mass energy region in terms 
of 5*, 1^0, and D2 partial waves. The D2 wave is dominated by the /2(1270) while the S wave contribution includes 
at least one additional resonance {fo{Y)) besides the /o(980), which could be the /o(1370) or the /o(1500). The fit 
includes the /o(980), another scalar, /o(^), /2(1270), and /2(1525). Known parameters are used for the /2(1270) 
and the /2(1525) in the fit. ro2 is the ratio of helicity to helicity 2 of the /2(1270). 
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The fit results are reported in Table V. The fit with the fo{Y) included is strongly favored. BES studied J/i/j ^ jTrn 
using 58 million J/V^ events and found a scalar with m = 1466 ± 6 ± 20 MeV/c^ and F = 108t\t ± 25 MeV/c^ [35], 
in good agreement with Belle's result. 
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FIG. 6: Partial wave analysis fit of the 27 to differential cross section in the low CM energy (W) region in terms of 5*, 

Do, and D2 partial waves by Belle [34]. The contributions of the components are also shown. 



TABLE V: Fit results for 27 production of tt^tt^ in the low mass energy region by Belle [34]. 



Parameter 


Nominal 


ro2 = 


No/o(y) 


Units 


mass(/o(980)) 


982.2 zb 1.0 


980.2 ± 1.0 


983.7l};o 


MeV/c^ 


r^^(/o(980)) 


285.5t}«3:? 


297.7li^;? 


370.5t?°;^ 


eV 


mass(/o(y)) 


1469.7 ± 4.7 


1466.8 =b 0.6 




MeV/c^ 


r(/o(K)) 


89.7l^-i 


499 4+18.4 




MeV 


r^^i3(/o(y) ^ 7r°^°) 


-'--'-•^-4.0 


6780.2t^^^;^ 


(fixed) 


eV 




^Q+0.24 
O.Ub>_Q 29 


(fixed) 




% 


B(/2(1270 ^ 77) 


1.57 ±0.01 


1 62+°-°2 
-'-•'^^-0.01 




xlO"^ 


X\ndf) 


1010 (615) 


1206 (617) 


1253 (619) 





V. BABAR ISR PHYSICS 

Because of the very high luminosity at B factories, hadron spectroscopy has also benefited greatly from studies 
using ISR to reduce the CM energy below the T(46') to study e+e~ hadrons from 1 < < 5 GeV. The BaBar 
collaboration have used 232 fb~^ at the T(46') to study e+e~ K~^K~7r^ and KsK~ 7r~^ using this technique. They 
require that the ISR 7 be detected, which forces the hadrons to be within the fiducial volume of the detector, and 
fully reconstruct the hadronic final state. The Dalitz plots are dominated by K^K* production. From a Dalitz plot 
analysis, they have separated the isoscalar and isovector components and measured their cross sections, as shown in 
Fig. 7. The channels are dominated by resonances that are consistent with the 0(1680) and p(1450) [36]. 

BaBar has also studied e~^e~ 7r^7r~7r^7r^ using ISR events [37]. Channels with four pions dominate the cross 
sections in the 1 < Ecm < 2 GeV region, and are very important for determinations of the anomalous magnetic 
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E_(G.V) E.„(0-V) 

FIG. 7: The cross section for (a) isoscalar and (b) isovector components from a Dalitz plot analysis of ISR e^e~ KKn by 
BaBar [36]. 

moment, a^, and the fine structure constant evaluated at the Z-pole, a{m'^). The cross section, shown in Figs. 8 a 
and b, is consistent with SND at low energy and is a huge improvement about 1.4 GeV, as shown in Fig. 8 b. The 
preliminary precision is 8%, and it is hoped that it will reach 5% over the peak region, which will help improve the 
precision of a^. This method will be used to improve the precision of R values at low energy, which Marco Verzocchi 
pointed out at this conference could be a bottleneck to future tests of Electroweak physics. 




* loou i2txi lAfiO me zmt 

E.. [GeVl E,, [MeVl 



FIG. 8: Cross section for ISR e^e tt^tt tv^tv^ [37]. (a) Full energy range, (b) Low energy range with results from previous 
experiments. Results are preliminary. 



VI. BEPCII AND BESIII 

Finally I report briefiy the status of BEPCII and BESIII. BEPCII is a two-ring e+e~ collider that will run in 
the tau-charm energy region {Ecm = 2.0 — 4.2 GeV, but possibly as high as 4.6 GeV) with a design luminosity of 
1 X 10^^ cm~^s~^ at a beam energy of 1.89 GeV, an improvement of a factor of 100 in luminosity with respect to the 
BEPC. This is accomplished mainly by using multi-bunches and micro-beta. 

The BESIII detector consists of a beryllium beam pipe, a helium-based small-celled drift chamber, Time-Of-Flight 
counters (TOF) for particle identification, a CsI(Tl) crystal calorimeter, a super-conducting solenoidal magnet with 
a field of 1 Tesla, and a muon identifier using the magnet yoke interleaved with Resistive Plate Counters. Fig. 9 
shows the schematic view of the BESIII detector, including both the barrel and end cap portions. 
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FIG. 9: Schematic view of the BESIII detector. 




FIG. 10: BESIII detector at the IP in June 2008. Shown are a superconducting quadrupole, the two beam lines, and the 
BESIII detector with the magnet iron open and end caps exposed. 



The detector moved to the IP in the spring of this year and is shown in Fig. 10 at its final location in June 2008 
with all beam magnets and vacuum pipes in place. Commissioning of the detector and collider together began in 
July, and the first hadronic event was obtained on July 19, 2008. Currently data at the '0(25') is being taken for 
calibration purposes. 

Clearly BESIII with a luminosity of 1 x 10^^ cm~^ s~^ will contribute greatly to precision fiavor physics; Vcd and 
Vcs will be measured with a statistical accuracy of better than 1.0%. D^D'^ mixing will be studied, and CP violation 
will be searched for. Huge J/il) and ip{2S) samples will be obtained. The rjc^ XcJi and he can be studied with high 
statistics. The high statistics will allow searches for physics beyond the standard model. The future is very bright. 
More detail on BEPCII and BESIII may be found in Ref. [38]. 



Acknowledgments 



The author wishes to thank the CLEOc, Belle, and BaBar collaborations for their help in preparing this paper 
and especially wants to thank Ryan Mitchell, Selina Li, Chengping Shen,and Xiaoyan Shen for their help. 



9 



34th International Conference on High Energy Physics, Philadelphia, 2008 
Work supported by Department of Energy contract DE-FG02-04ER41291. 



[1] E. Eichten et al, Rev. Mod. Phys. 80, 1161 (2008). 

[2] H. Mendes et a/., CLEO collab., Phys. Rev. D 78, 011102 (2008). 

[3] PDG: C. Amsler et a/., Phys. Lett. B 667, 1 (2008). 

[4] Dudek et a/., Phys Rev. D 73, 07450 (2006). 

[5] R. E. Mitchell et a/., CLEO collab., arXiv:0805.0252 (2008). 

[6] S. Uehara et al, Belle Collab., Eur. Phys. Jour. C 53, 1 (2008). 

[7] H. Nakazawa, Belle Collab., Talk at Photon 2007, Paris 9 - 13 July, 2007. 

[8] C. Davies, Talk at 34th International Conference on High Energy Physics, Philadelphia, Pennsylvania, (2008). 

[9] M. Andreotti et al, E835 Collab., Phys. Rev. D 72, 032001 (2005). 

[10] J. L. Rosner et a/., CLEO Collab., Phys Rev. Lett. 95, 102003 (2005); P. Rubin et al, CLEO Collab., Phys. Rev. D 72, 
092004 (2005). 

[11] S. Dobbs et a/., CLEO Collab., arXiv:0805.4599 (2008). 

[12] V. A. Novikov, Phys. Rep. 41, 1 (1978). 

[13] S. K. M. Ecklund et a/., CLEO Collab., arXiv:0803.2869 (2008). 

[14] L. Landau, Phys. Abstracts A 52, 125 (1949); C. N. Yang, Phys. Rev. 77, 242 (1950). 

[15] R. Barbieri et a/., Phys. Rev. D 37, 3210 (1988); Nucl. Phys B 192, 61, (1981). 

[16] M. Ablikim et al, BES Collab., Phys. Rev. Lett. 97, 121801 (2006). 

[17] M. Ablikim et a/., BES Collab., Phys. Lett. B 641, 145 (2006). 

[18] M. Ablikim et al, BES Collab., Phys. Rev. D 76, 122002 (2007). 

[19] M. Ablikim et al, BES Collab., Phys. Lett. B 659, 74 (2008). 

[20] J. Z. Bai et a/., BES Collab., Phys. Lett. B 605, 63 (2005). 

[21] N. E. Adam et al, CLEO Collab., Phys. Rev. Lett. 96, 082004 (2006); T. E. Coans et al, CLEO Collab., Phys. Rev. 
Lett. 96, 182002 (2006); G. S. Huang et a/., CLEO Collab, Phys. Rev. Lett 96, 032003 (2006); G. S. Adams et a/., CLEO 
collab., Phys. Rev. D 73, 012002 (2006); D. Cronin-Hennessy et al, CLEO Collab., Phys. Rev. D 74, 012005 (2006); R. 
A. Briere et a/., CLEO Collab., Phys. Rev. D 74, 031106 (2006). 

[22] M. Ablikim et a/., BES Collab., arXiv0807.0494 (2008). 

[23] B. Aubert et al, BaBar Collab., Phys. Rev. D 74, 091103 (2006). 

[24] B. Aubert et al, BaBar Collab., Phys. Rev. Lett. 95, 142001 (2005). 

[25] M. Ablikim et a/., BES Collab., Phys. Rev. Lett. 100, 102003 (2008). 

[26] L Adachi et aL, Belle Collab., arXiv:0808.0006 [hep-ex] (2008). 

[27] F. Mane et a/., DM1 Collab., Phys Lett. B 112, 178 (1982). 

[28] B. Aubert et aL, BaBar Collab., Phys. Rev. D 77, 092002 (2008). 

[29] G. J. Ding and M. L. Yan, Phys. Lett. B 650, 390 (2007). 

[30] G. J. Ding and M. L. Yan, Phys. Lett. B 657, 49 (2007). 

[31] Z. G. Wang, Nucl. Phys. A 791, 106 (2007). 

[32] E. Klempt and A. Zaitsev, Phys. Rep. 454, 1 (2007). 

[33] T. Mori, Belle Collab., arXiv:0704.3538 (2007). 

[34] S. Uehara, Belle Collab., accepted by Phys. Rev. D, arXiv:0805.3387 (2008). 

[35] M. Ablikim et aL, BES Collab., Phys. Lett. B 642, 441 (2006). 

[36] B. Aubert et aL, BaBar Collab., Phys. Rev. D 77, 092002 (2008). 

[37] S. Li, Talk at 34th International Conference on High Energy Physics, Philadelphia, Pennsylvania, arXiv 0810.0796 (2008). 

[38] F. A. Harris, BES Collab., Talk at the 10th International Workshop on Meson Production, Properties and Interaction 
(Meson08), Cracow, Poland, 6-10 June, 2008, arXiv0808.3163 (2008). 



10 



